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The first order phase transition proceeds via nucleation and growth of 
true vacuum bubbles. When charged particles collide with the bubble they could 
radiate electromagnetic wave. We show that, due to an energy loss of the particles 
by the radiation, the damping pressure acting on the bubble wall depends on the 
velocity of the wall even in a thermal equilibrium state. 
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There have been many studies on cosmological roles of first-order phase tran- 
sitions which proceed via nucleation and growth of vacuum bubbles. Studying 
particle scatterings at a moving bubble wall is important to know the bubble 
dynamics in a hot plasma. For example, to calculate the velocity of electro-weak 
bubbles |l|, 0] and the CP violating charge transport rate by the wall available for 
baryogenesis, one should know the reaction force acting on the wall by the 
scattered particle(e.g'. quarks, gauge bosons). 

When charged particles collide with a bubble wall, they could emit electro- 
manetic radiation. In this paper we study the effect of energy being carried 
away by the radiation on the pressure acting on the bubble wall. When particles 
enter a true vacuum bubble through a wall (say, the electroweak bubble), they 
could acquire a mass and are decelerated. For example, a fermion ijj may get a 
mass through the well known Yukawa term gilxpip = mipil), where is a Higgs 
field. At this time if the particle is charged electomagnetically, it can radiate 
electromanetic wave due to the deceleration. Let us calculate the pressure that 
scattered particles exert on the bubble wall. For simplicity we assume a linear 
profile for the bubble wall. (See Fig.l.) and choose a rest frame of the bubble 
wall. Then,m(x) = rriQx/d {0 < x < d) is the position dependent rest mass of 
the particle at the wall. We also assume that the mean free path of the collision 
is much shorter than the wall width d, hence WKB approximation is good. The 
radiation power of an accelerated particle is given by the relativistic version of 
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Lamor's formula: 

/ -\ 2 

dErad _ a I dk 
dt m? I dt 

where a = 2e^/3c^ ~ 0.0611 in natural units h = c = k = 1 and k = (k^, ky, k^) 
is a 3-momentum of the particle. Assuming that the wall is planar and parallel 
to the y — z plane we can treat the bubble as 1-dimensional one along the x-axis. 
Energy, momentum and mass of the particle satisfy an usual relation 

^2 = m2(x) + |A?p(x). (2) 

Let us denote k^ as k from now on. Differentiating the above equation with time 
t and using k = Edx/dt, we get a force 

dk drn^ 1 

~dt ^ ~~d^2E' ^ ' 

which is a starting point of the pressure calculation [Q. However, if we also con- 
sider the energy carried away by the radiation Erad, then the total energy con- 
served is Etot = Erad+E and the force and , hence, the pressure should be changed. 
From dEtot/dt = we obtain 

k dk mdm a f dk\^ 



which has a solution 



dk —m 
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dt 2aE 



_^ I AaEk dm 
~ V m4 dt 



(5) 



Up to 0{a) we can expand the root term and get 

dk dmm 2a I dm\ 
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dt dx E kE \dx J 



(6) 



Then the total pressure by the colhsion of the particles in the plasma is given 
byi 



P= dx 



d^k 



(7) 



(27r)3 

where f{E) = {Exp{(3E) ± 1)"^ is a distribution function of fermions and bosons 
respectively. First, let us briefly review the well-known results without the radi- 
ation damping( a = 0). When the mean velocity of the plasma fluid V (or minus 
of the bubble wall velocity relative to the fluid ) is zero, the first term of the 
Eq.® contributes 

^ r°° dm^(x) , r d^k 1 1 

Pi = ; dx 



dx J (27r)3 2Ee'3^±l 
= F(m,r)-F(0,T), (8) 

where F{(j), T) is the T dependent part of free energy at a temperature T = (3~^. 
When V ^ Q the distribution function is changed to 

/[7(i?-\/fc)] = (e^^(^-^'=)±l)"\ (9) 

where 7 is the Lorentz gamma factor. However, using the fact that the phase 
factor d^k/E is Lorentz invariant and changing the integration variable to k' = 
7(A; — VE) and defining E' = 'y{E — Vk) one can find that the V dependency of 
Pi disappears 0. So one have needed to consider a non-equilibrium deviation of 
/ to calculate the velocity of the wall 0. 

Now let us consider the effect of the radiation, the second term of Eq.(B]). 
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When V = the term contributes 



oo 



which vanishes, because the integrand is an odd function of k. When V ^ one 
can easily check that due to the 1/k term the V dependency survives even under 
the change of the integration variable. So in this case 



°° ^ (im(x)\^ f d^k 1 



Po = 2a — ^ dx 



oo 



oo 



= 2a I ( ^^^1 dx h{x) (11) 



dx ) J (27r)3Efc(e/57(s-^^)±l) 
77,(2 
00 \ dx J 

Let us discuss overall features of this integration qualitatively. To investgate 
a singular behavior of the integrand in Eq. ([Tl|) near A; = we first consider k 
integration(along x-axis) 

/ . r iR,, = r iiJpMdk + ,(0) r Uk, (12, 

J — oc K J — 00 rC — U J — 00 rC 

where g{k) = {E{Exp[j3'~f{E — Vk)] ± l})^"*^. The second term of the above 
equation is zero. For the first term, since the contribution from the integrand 
with ~ dominates, one can approximate the term as 

r ^dkc^g{ko)~g{-ko), (13) 
J-ko dk 

where [ko, —ko] is a small momentum interval in which the integrand contributes 
significantly. If \^ = two terms cancel each other. When y <^ 1 one can expand 
them up to 0{V) and expect that they are linearly dependent on V. On the other 



hand, as ^ 1 (7 — 00) they are exponentially suppressed. The integration 
over ky and kz are just sums of / with a shifted value of E, hence the overall 
behaviour would not be changed after 3- dimensional calculation. ( During the 
numerical study it is useful to change the measure from dkydk^ to 2nEdE) The 
numerical integration of Eq.(0) confirms this behaviour. (See Fig. 2) 

To be more concrete, let us obtain an approximate value of the integration 
when V 1 for fermions. In this case we can expand /[7(-E — Vk)] ~ f{E) — 
Vf3kf{E)[f{E) - 1] = f{E) + V(3kp{E)Exp{-f3E). The integration with the 
first term gives zero and the second term contributes 

, because 

to lowest oder in (m/T)^ (See Ref.p). Therefore, for the wall described in Fig.l, 
the pressure by the radiation is 

P. = (M^^, (16) 

which well fits with the numerical result shown in Fig. 2 when V is small. This 
pressure is proportional to the wall velocity up to the moderately relativistic 
case and exists even when the system is in a thermal equilibrium. Since P2 is 
proportional to the square of the charge, an antiparticle contributes equally as a 
particle. 



In summary we consider the radiation damping of the particles colhding with 
bubbles and calculate the velocity dependent pressure on the bubble wall. 

Authors are thankful to Dr. Myongtak Choi for useful discussions. This 
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FIG.l. 

The mass of the particle m{x) in the wall rest frame. 



FIG.2. 

The pressure by the radiation damping of fermions colliding with the linear bubble 
wall (Eq.(pJ])). Here 1/d = 1 = mo = T for simplicity. 
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Fig.l. 
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